Divergent Precursors of the Mott-Hubbard Transition at the Two-Particle Level 
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Identifying the fingerprints of the Mott-Ifubbard metal-insulator transition may be quite elusive 
in correlated metallic systems if the analysis is limited to the single particle level. However, our 
dynamical mean-field calculations demonstrate that the situation changes completely if the frequency 
dependence of the two-particle vertex functions is considered: The first non-perturbative precursors 
of the Mott physics are unambiguously identified well inside the metallic regime by the divergence 
of the local Bethe-Salpeter equation in the charge channel. At low temperatures this occurs in 
the region where incoherent high-energy features emerge in the spectral function, while at high 
temperatures it is traceable up to the atomic-limit. 
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Introduction. Among all fascinating phenomena charac- 
terizing the physics of correlated electronic systems, one 
of the most important is undoubtedly the Mott-Hubbard 
metal-to- insulator transition (MIT)|l|. Here, the onset of 
an insulating state is a direct consequence of the strong 
Coulomb repulsion, rather than of the underlying elec- 
tronic bandstructure. Mott MITs have been indeed iden- 
tified in several correlated materials Q, especially in the 
class of transition metal oxides and heavy fermions. The 
interest in the Mott MIT is not limited however to the 
transition "per se" , but it is also for the correlated (bad) 
metallic regime in its proximity. In fact, this region of the 
phase-diagram often displays a rich collection of intrigu- 
ing or exotic phases, that many relate also to the physics 
of the high-temperature superconducting cuprates. 

An exact theoretical description of the Mott MIT rep- 
resents a considerable challenge due to its intrinsically 
non-perturbative nature in terms of the electronic in- 
teraction. However, a significant progress was achieved 
with the invention of the dynamical mean field theory 
(DMFT)[i]. By an accurate treatment of local quan- 
tum correlations, DMFT has allowed for the first non- 
perturbative analysis of the Mott-Hubbard MIT in the 
Hubbard model[6[, and, in combination with ab-initio 
methods f^, also for the interpretation and the prediction 
of experimental spectroscopic results for strongly corre- 
lated materials, such as, e.g., the paramagnetic phases of 
V2O3. Theoretically, several "hallmarks" of the onset of 
the Mott insulating phase can be unambiguously identi- 
fied in DMFT: At the one-particle level, a divergence of 
the local electronic self-energy in the zero-frequency limit 
is observed, reflecting the opening of the Mott spectral 
gap, while, at the two-particle level, the local spin sus- 
ceptibility {xs{^^ = 0)) diverges at T = 0, due to the onset 
of long-living local magnetic moments in the Mott phase. 
Description of the problem. While the characterization 
of the MIT itself is quite clear, at least on a DMFT level, 



the physics of the correlated metal regime in the vicin- 
ity of the MIT is far from being trivial and presents 
several anomalies. We recall here: the occurrence of 
kinks of purely electronic origin Q in the angular re- 
solved one-particle spectral functions or in the electronic 
specific heat, the formation of large instantaneous mag- 
netic moments, screened by the metallic dynamics^, the 
abrupt change of the out-of-equilibrium behavior after a 
quench of the electronic interaction [l^l, and the changes 
in the energy-balance between the paramagnetic and the 
low-temperature (antiferromagnetically) ordered phase, 
which also affect the restricted optical sum-rules [11 1. 
Partly motivated by these observations, many DMFT 
calculations have been aiming at general characterization 
of this regime, e.g., by studying the phase-diagram of the 
half-filled Hubbard model. However, no trace of other 
phase transitions has been found beyond the MIT itself 
and the (essentially) mean-field antiferromagnetically or- 
dered phase, which is not of interest here. Hence, one 
of the main outcomes of the preceding DMFT analyses, 
mostly focusing on the evolution of one-particle spectral 
properties (and, to less extent, on susceptibilities[12]), 
has been the definition of the "borders" of the so-called 
crossover regions at higher T than those where the MIT 
can be observed. The shape of these crossover regions 
has been analyzed in many different waysdIiMI- We 
note here, that the (different) criteria used for defining 
crossover regimes imply a certain degree of arbitrariness. 
Furthermore, the crossover region is located at much 
higher Ts than those, where some of the abovementioned 
anomalies are observed. 

In this paper, going beyond the standard, typically 
one-particle, DMFT analyses, we present a completely 
unambiguous criterion to distinguish the "weakly" and 
the "strongly" -correlated regions in the phase-diagram. 
By studying the frequency structure of the two-particle 
local vertex functions of DMFT, we observe the diver- 
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gence of the local Bethe-Salpeter equation in the charge 
channel. This divergence defines a regime remarkably 
different, also in the shape, from the crossover region, 
where non-perturbative precursor effects of the MIT be- 
come active, even well inside the low-temperature metal- 
lic phase. The precise definition of such a regime rep- 
resents a much better playground for a general interpre- 
tation of the anomalous physics emerging as a precursor 
of the MIT. Furthermore, our analysis, showing the oc- 
currence of peculiar divergent features in some of the 
two-particle local vertex functions of DMFT is also ex- 
pected to have a significant impact on future calculations 
for strongly correlated electron systems, because the two- 
particle local vertex functions represent a crucial ingredi- 
ent for both (i) the calculation of dynamical momentum- 
dependent susceptibilities in DMF T [3. [ItI . [18{ , as weh 
as (ii) the diagrammatic extensions 19|,|20| of the DMFT, 



aiming at the inclusion of non-local spatial correlations. 
DMFT results at the two particle-level. We consider the 
Hubbard model on a square lattice in the paramagnetic 
phase at half-filling, that is one of the most basic realiza- 
tions of the MIT in DMFT. The corresponding Hamilto- 
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where t is the hopping amplitude between nearest- 
neighbors, U is the local Coulomb interaction, and 
cl^{cia) creates (annihilates) an electron with spin a =t 
, J, at site i; = cl^Cia. Hereafter, all energy scales will 
be given in units oi D ~ At ^ 1, i.e., of the half of the 
standard deviation of the non-interacting D0Si21]. 

Differently from previous studies, we will focus on the 
analysis of the two-particle local vertex functions com- 
puted with DMFT. By using an Hirsch-Fye quantum 
Monte Carlo impurity solver [J|, whose accuracy has been 
also tested in selected cases with exact-diagonalization 
DMFT calculations, we have first computed the general- 
ized local susceptibility ^'^'^ (w). This is defined, follow- 



ing the notation of Ref. [22|, as: 



dTidT2dT3 e-*'^^ie*('^+'^)^^e-*('''+'^)^^ (2) 
X (T,4(Ti)c,(r2)4,(T3)c,-(0)) 
- (r.4(ri)c.(r2))(T.4,(r3)c.,(0))' 

where T^- is the (imaginary) time ordering operator, and 
v' and uj denote the two fermionic and the bosonic 
Matsubara frequencies, respectively. Then, the Bethe- 
Salpeter equation in the charge channel (defined as 
Xc"^ (i^) — x'f'f {^) + X'fi {^) has been considered for 
w = 0, allowing to determine the corresponding irre- 
ducible vertex 1221 



[xr {io^or'-ix^o'' (^ = 0)]- 
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where the last term is defined through the convolu- 
tion of two DMFT Green's functions, as Xo'^ {'^) = 
-T-^G{v)G{iy+uj)Si,^' . The vertex F^'"' can be viewed as 
the two-particle counter-part of the electronic self-energy 
and, for a half-filled system, it is a purely real function. 
Our numerical results are reported in Fig. [T] for four dif- 
ferent values of the electronic interaction ?7 at a fixed 
temperature of T = 0.1. Starting to examine the first 
panel, corresponding to the smallest value of C/ = 1.2, 
one observes two main diagonal structures in the Mat- 
subara frequency space: These structures are easily in- 
terpretable as originated by reducible ladder-processes in 
the (transverse) particle- hole {i' = iy') and in the particle- 
particle {ly = —ly') channels respectively [22,]. Following 
the behavior of the local spin susceptibility in the Mott 
phase, the main diagonal structure will diverge exactly 
at the MIT (Uun ~ 3) in the T = limit. 

In contrast to these standard properties of V^" , visible 
in the first panel, the analysis of the other three panels 
of Fig. [T] shows the emergence of a low-frequency sin- 
gular behavior of the vertex functions for a value of U 
much smaller than that of the MIT: Already aXU = 1.27 
(second panel), one observes a strong enhancement of 
the vertex function at the lowest Matsubara frequencies 
(note the change in the intensity-scale). This is visible 
as an emergent "butterfly" -shaped structure, where the 
intense red-blue color coding indicates alternating signs 
in the {vjv') space. Remarkably, such a low-energy struc- 
ture becomes predominant over the other ones along the 
diagonals. That a true divergence takes place is sug- 
gested by the third panel [U = 1.28), where the intensity 
of the "butterfly" structure is equally strong but the signs 
are now inverted as indicated by the colors. This is also 
shown more quantitatively by a selected cut of FJ^" in 
frequency space, reported in the corresponding panels of 
the second row of Fig. 1. Note that the inversion of the 
signs cannot be captured by perturbation theory calcu- 
lations (green circles), marking the non-perturbative na- 
ture of the result. The rigorous proof of the divergence 
is provided by the evolution of the matrix Xc^ j which is 
positive deflnite at weak-coupling, while one of its eigen- 
values (see legends and insets in the bottom row of Fig. [1]) 
becomes negative crossing 0. Finally, by further increas- 
ing J7, the low-energy structure weakens, indicating that 
at fixed T = 0.1 this vertex divergence is taking place 
only for a speciflc value of the Hubbard interaction, i.e., 
for U ~ 1.275. This flnding naturally leads to the crucial 
question of the temperature dependence of the results: 
Does such divergence occur for all temperatures, and if 
yes, is the temperature dependence of U signiflcant? As 
one can immediately understand from FiglS] the answer 
to both questions is positive 23]: By repeating the anal- 
ysis of FigUat different temperatures, we could identify 
the loci (T, red dots in Figl^]) in the phase-diagram, 
where the low-frequency divergence of FJ^"^ occurs. This 
deflnes a curve T(U) with a quite peculiar shape, where 
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FIG. 1. (Color Online) Upper row: Evolution of the frequency dependent two-particle vertex function, irreducible in the 
charge channel, (F"" ) for increasing U. The data have been obtained by DMFT at zero external frequency (a; = 0) and fixed 
temperature (T = 0.1); lower row: linear snapshot of the same Fc along the path marked by the dashed-line in the first panel 
of the upper row, i.e., as a function of i/ = j{2n + 1) for n' = {u' — ^), compared to perturbation theory (FT) results. In 

the legends/insets the closest-to-zero eigenvalue (A) of Xc'" /Xo " is reported for each U. 



three regions can be distinguished: (I) at very high T, 
the behavior is almost perfectly linear T cx [/ ; (II) in the 
low T limit the curve strongly bends, extrapolating for 
T ^ at U{0) - 1.5 < [/mit ^ 3; (III) at intermediate 
T the curve interpolates between these two regimes, with 
a "re-entrance" clearly affected by the presence of the 
MIT at larger U (blue squares in Fig. [2]). We note that 
by increasing U much further than the T{U) curve, one 
eventually observes a divergence also of the local Bethe- 
Salpeter in the particle-particle channel (orange points in 
Fig. [2|), while for all values of T, U considered, no similar 
divergence is found in the spin channel. 

Interpretation of the results. In contrast to the case of 
the main diagonal structures of the vertex functions, the 
interpretation of the low-frequency divergences of TJ^'^ 
is not directly related to the MIT. However, even if at 
low T the divergences take place in the metallic region 
of the phase-diagram, the re-entrance shape of the T{U) 
curve is indeed remarkably affected by the position of 
the MIT: The most natural interpretation is, hence, that 
the shaded area in the phase diagram defines the region 
where the precursor effects of the MIT physics preclude 
the perturbative description and become a crucial ingre- 
dient in determining the properties of the system. This 
interpretation is evidently supported by the fact that the 
signs of the two-particle vertex functions are correctly 
predicted in perturbation theory only up to the left-hand 
side of the T(U) curve. More generally, the T{U) curve 



can be identified as the limit of the region of applica- 
bility of schemes based on the Baym-Kadanoff 24| func- 



tional ^[G], since |^ — Fc can no longer be defined 
on that line. At the same time, the low-frequency singu- 
larities of the vertex may render problematic the numeri- 
cal evaluation of the Bethe-Salpeter equation to compute 
momentum-dependent DMFT response functions in spe- 
cific regions of the phase-diagrams, suggesting the use of 
alternative procedures [25|. 

As the singularity of Fc (and later on of Tpp) is not 
associated to simultaneous divergences in the other chan- 
nels, the application of the local parquet equations [iii [26l | 
allows to identify the ultimate root of these divergences in 
the fully two-particle irreducible diagrams. Hence, this is 
an "intrinsic" divergence, deeply rooted in the diagram- 
matics and not generated by ladder scattering processes 
in any channel. From a more physical point of view, 
the fact that the only irreducible vertex F displaying 
no-singularities at low-frequencies is the spin one might 
also indicate the emergent role played by preformed lo- 
cal magnetic moments as MIT precursors, even in regions 
where the metallic screening is rather effective. 

We can go, however, also beyond these general consid- 
erations and analyze the three regimes of the T{U) curve 
in detail, also discussing the relation with the emergence 
of some of the anomalous properties of the physics in the 
vicinity of the MIT. The analysis of the high-T linear 
regime [(I) in Fig. [5] oiT{U) is probably the easiest: Here 
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FIG. 2. (Color Online) Instability lines of the irreducible ver- 
tices in the charge (Fc red circles) and in the particle-particle 
channels {Tpp orange diamonds) reported in the DMFT phase 
diagram of the half-filled Hubbard model (the data of the 
MIT, blue solid line, are taken from Ref [il,[i3). The red 
dashed line indicate the corresponding instability condition 
(T = -^U) estimated from the atomic limit. Inset: zoom 
on the low-T region, where also different estimations (dashed 
light-blue[l3|, dashed bluep^]) of the crossover region are in- 
dicated. 



U,T ^ D, and hence a connection with the atomic limit 



0) can be done: Using analytic expressions [22, l28| 
for the reducible two-particle vertex functions as an in- 
put for Eq. dS]), we find that the low- frequency divergence 
of T'^'^' occurs at T/U = ^ and that the eigenvector 
associated to the vanishing eigenvalue of Xc'^ has the 
particularly simple form: -^{S^f^rp-f — K{-ttT))- As it is 
clear from the comparison with the red dashed line in 
Fig. [2] this proportionality exactly matches the high-T 
linear behavior of our TiU) curve. Crossing this curve 
in its high-T linear regime, which extends indeed over a 
large portion of the phase-diagram, corresponds to enter- 
ing a region where the thermal occupation of the high- 
energy doubly-occupied/empty states becomes negligi- 
ble, letting the physics be dominated by the local mo- 
ments. The connection with the local moment physics 
also holds for the low-T region (II), though via a differ- 
ent mechanism: For T — >■ 0, the relevant energy scales are 
the kinetic D) and the potential {U) energy, whose 
competition is regulated by quantum fluctuations. In 
this case, obviously, only numerical results are available: 
We observe that the extrapolated value of f/(0) 1.5 
falls in the same region (gray arrow in inset of Fig. [2]), 
where DMFT(NRG)[2'7l see a first clear separation of 
the Hubbard sub-bands from the central quasi-particle 
peak in the spectral function A{uj). We recall here that 
the formation of well-defined minima in A(w) between 
the central quasi-particle peak and the Hubbard sub- 
bands, is directly connected with the anomalous phe- 
nomenon of the appearance of kinks in the electronic 



self-energy and specific heat[8|. At the same time, more 
recent^DMFT(DMRG)t2i| data rather indicate that for 
U >U{Q) ^ 1.5 two sharp peak-features emerge at the 
inner edges of the Hubbard sub-bands, which, however, 
would be already visible aiU >1. 

Looking for a more analytical description of this sce- 
nario, we can consider the DMFT solution of the much 
simpler Falicov-Kimball (FK) model[30]: Here one can 
exactly show that F^"^ indeed diverges before the MIT is 
reached (precisely at: U^^ — -^U^r^). However, for the 
FK results, a direct relation with the formation of the two 
minima in Aiuj) cannot be completely identified, as the 
renormalization of the central peak is not captured in this 
scheme ^3]. Finally, an interesting observation can be 
made about the most complicate intermediate T regime 
(III), where all energy scales (D, U, T) are competing: 
Recent out-of-equilibrium calculations for the Hubbard 
model have shown [l3|, that after a quench of the inter- 
action (i.e., from U = to J7 > 0), the system's relax- 
ation occurs in two different (non-thermal) ways. The 
changeover between these two regimes, however, appears 
for a given set of parameters U ~ 1.65 and T^s ~ 0.4, in 
close proximity of the vertex instability-line in our phase 
diagram. 

Conclusions and outlook. Our DMFT calculations have 
shown how the emergent (non-perturbative) precursor ef- 
fects of the MIT determine a low-frequency divergence 
of the local Bethe-Salpeter equation in the charge chan- 
nel. This allows for an unambiguous identification of the 
regime, where perturbation and Baym-Kadanoff func- 
tional theory break down, and where, at the same time, 
several anomalous properties are observed or predicted 
for correlated metals. Taking properly into account the 
physics emerging from the singularity of the two-particle 
vertex functions will represent one of the main challenges 
for future improvements of the theoretical many-body 
treatments at the precision level required by the increas- 
ingly higher experimental standards. 
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